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ABSTRACT 


P 


In  a  previous  paper,  Jewell  and  Sundt  showed  how  to  approximate  the 
distribution  of  total  losses  from  a  large,  fixed,  heterogeneous  port¬ 
folio,  using  a  recursive  algorithm  developed  by  Panjer  for  the  dis¬ 
tribution  of  a  random  sum  of  random  variables  (a  single  casualty 
contract).  This  paper  extends  the  approximation  procedure  to  large, 
dynamic  heterogeneous  portfolios,  in  order  to  model  either  a  portfolio 
of  correlated  casualty  contracts,  or  a  future  portfolio,  whose  composi¬ 
tion  is  not  known  with  certainty. 


APPROXIMATING  THE  DISTRIBUTION  OF  A  DYNAMIC  RISK  PORTFOLIO 


by 

William  S .  Jewell 


0.  INTRODUCTION 

The  problem  of  finding  the  distribution  of  y  -  x^  +  +  . . .  +  Xjj  , 

where  the  (x^)  are  a  fixed  and  large  set  of  independent,  nonidentically 
distributed,  integer-valued  random  variables  was  considered  in  Jewell  and 
Sundt  (1981)  (hereinafter  referred  to  as  JS).  Although,  in  theory,  the 
discrete  density  of  y  is  just  the  N-fold  convolution  of  the  Individual 
densities,  this  computation  is  very  time-consuming,  and  various  forms 
of  approximation  must  be  used;  moreover,  in  many  risk  applications,  the 
use  of  a  normal  approximation  gives  very  bad  results,  even  for  large  N  , 
because  of  the  skewness  and  long  tails  of  the  density.  However,  if  the 

probability  p^  -  Pr  {x^  -  0}  is  significant  for  most  i  -  1,2 . N  , 

it  turns  out  that  a  very  good  approximation  can  be  obtained  using  newly- 
developed  procedures  for  the  related  problem  of  calculating  the  distribu¬ 
tion  of  the  sum  of  a  random  number  of  independent  and  identically  dis¬ 
tributed  random  variables. 

In  many  risk  applications,  especially  in  insurance  and  investment 
management,  there  are  an  ever-changing  number  of  risks  of  different  types 
and  it  is  of  interest  to  predict  the  distribution  of  a  portfolio  whose 
future  composition  is  not  known  with  certainty.  This  paper  develops  a 
general  model  for  this  situation,  and  shows  how  the  approximation  proce¬ 
dure  described  in  JS  can  be  extended. 


t-V. 


1.  THE  DYNAMIC  PORTFOLIO  MODEL 

Let  i  *  1,2,  . N  Index  a  number  of  different  risk  classes 
(Insurance  policies  or  types  of  investment)  in  a  given  portfolio,  and  let 
e  [0,1,  ...]  be  the  random  number  of  independent  risks  of  type  1  , 
giving  a  grand  total  number  of  risks  in  the  portfolio 

N 

(1-1)  iL  *  I  n  . 

1  i-1  1 

Risks  of  type  i  are  similar ,  in  the  sense  that,  if  is  the 

random  monetary  gamble  from  the  risk  of  type  1  ,  then  its  discrete 
density,  f°(x)  ,  is  the  same  for  all  j  ,  i.e.: 

(i  -  1,2 . N) 

(1.2)  Pr  {x,  -  x)  -  f?(x)  . 

J  (J  •  1,2 . n±) 

Ue  shall  only  consider  discrete  gambles,  with  the  common  range  of  the 
(x^j)  as  [0,1,2,  . R]  .  As  mentioned  above,  we  assume  for  the  moment 
that  the  (x^)  are  statistically  Independent  of  each  other  and  the  (n^)  , 
but  we  do  not  assume  that  the  (n^)  are  independent.  (But  see  Appendix  A.) 

The  total  monetary  gamble  for  all  risks  of  type  i  is  then  the  sum  of 
a  random  number  of  random  variables: 

{0  (n±  -  0) 

(i  -  1,2 . N) 

iu +  i12  +  ...  +  ilgi  >  0)  . 


and  the  grand  total  monetary  risk  is  then  the  fixed-term  random  sum: 


Note  that  the  (x^)  are  now  dependent  random  variables,  if  the  (n^) 
are. 

If  g(y)  and  iKn^n^  n^)  are  the  discrete  densities  of 

the  total  risk  sum  and  the  number  of  risks  of  each  type,  respectively, 
we  have  then  the  discrete  density  of  y  as: 


Pr  {y  -  y}  -  g(y)  -  \  \  ...  I  irfn^n^ 

nl  n2  °N 


(1.5) 


(y  “  0»1*2,  •••) 


which,  of  course,  is  a  lengthy  and  laborious  computation.  (In  JS,  the 
special  case  of  (n  )  deterministic  was  considered.) 


2.  INTERPRETATIONS 


Before  describing  a  method  of  approximating  (1.5),  we  give  some 
possible  practical  interpretations  of  the  model. 

In  insurance  applications,  the  simplest  interpretation  is  that  i 
refers  to  different,  distinguishable  types  of  insurance  policies  in  a 
given  portfolio;  for  instance,  similar  policies  in  personal  lines  of 
Insurance  could  refer  to  ordinary  life  insurance  policies  with  the 
same  face  values  issued  to  persons  of  the  same  age.  For  the  current 
year,  we  know  exactly  the  number  of  policies  of  type  i  and  hence, 
following  JS,  can  find  an  approximation  to  the  current  g(y)  .  How¬ 
ever,  an  approximation  to  (1.5)  would  be  necessary  to  predict  total 
portfolio  risk  for  next  year,  after  some  policies  are  withdrawn, 
some  policies  have  paid  out  benefits,  or  new  policies  have  been  added, 
and  still  others  have  shifted  type.  By  specifying  the  stochastic  law 
governing  this  "drop-add"  mechanism,  we  can  get  ...»  n^) 

for  next  year.  Possible  reasons  for  leaving  correlation  between  the 
(n^)  are  that  we  may  have  a  precise  idea  of  how  new  sales  are  dis¬ 
tributed  among  the  different  types  of  policies,  but  may  be  uncertain 
about  the  total  new  business;  or,  the  new  business  total  may  be 
accurately  estimated  but  the  distribution  may  be  uncertain;  or, 
there  may  be  an  uncertain  number  of  policies  which  are  shifting  type 
(as  in  aging  of  life  insurance  insured);  etc. 

A  second  insurance  interpretation  is  the  so-called  casualty  claim 
model,  in  which  multiple  claims  may  occur  on  a  policy  during  a  given 
exposure  year.  Here  i  Indexes  each  of  a  fixed  number  of  policies, 
f°(x)  is  the  individual  claim  ("severity")  density,  n^  is  the  random 


number  of  claims  ("frequency") ,  and  x^  is  now  the  total  monetary  claim 
on  the  single  policy  i  .  Of  course,  if  the  (n^)  were  independent, 

then  this  application  could  be  handled  by  making  £  ^(n^) 
the  basic  density  used  in  the  procedure  described  in  JS;  but  this  would 
require  prior  calculation  of  this  compound  law  (see  also  Section  8 
and  9,  below).  Moreover,  external  factors,  such  as  weather  and  economics, 
often  affect  the  number  of  claims  of  all  types  of  contracts  in  a  given 
portfolio  in  the  same  way,  thus  Introducing  correlation  and  the  need 
for  a  more  general  model. 

In  most  insurance  portfolios,  a  great  deal  of  effort  is  used  to 
assure  that  the  (x^)  are  statistically  independent  of  each  other. 

However,  there  remains  always  the  possibility  that  risks  of  the  same 
type  i  are  influenced  by  the  same  exogenous  factors.  In  Appendix  A, 
we  consider  the  case  when  risks  of  the  same  type  are  exchangeable  random 
variables,  which  leads  to  a  weak  form  of  dependence  on  the  (x^ )  . 

In  investment  portfolios,  it  is  unusual  to  have  Independent  risks 
of  the  same  type,  i.e.,  requiring  the  same  investment  level,  and  having 
the  same  outcome  distribution;  instead,  we  usually  have  a  different 
amount  of  money  invested  in  different  risks.  If  we  let  n^  be  the  level 
of  investment  in  type  i  and  x^  the  net  return  from  this  investment, 
then  (1.3)  holds  only  if  the  (x^)  are  perfectly  correlated,  or  what 
is  the  same,  if  (1.3)  is  replaced  by  x^  -  n^x^  .  Another  limitation 
on  Investment  modelling  is  that  it  is  usually  possible  to  have  negative  net 
returns,  which  is  discussed  in  Section  10.  It  should  be  remembered  also 
that  our  approximation  is  usually  successful  only  if  the  problem  is  modelled 
so  that  the  probability  of  zero  net  return  is  substantial;  i.e.,  all  "sure 
thing"  return  has  been  eliminated. 


Technological  risk  applications  are  based  upon  the  compound  law 
interpretation;  for  instance,  in  reliability  engineering,  n^  may  refer 
to  the  random  number  of  mechanical,  electrical,  or  thermal  shocks  of 
type  i  which  affect  a  given  piece  of  equipment;  in  fire  damage  analysis, 
is  the  random  number  of  fires  of  a  given  type  (size,  type  of  dwelling 
or  land  classification)  which  occur;  and  so  forth.  In  technology  appli¬ 
cations,  the  primary  modelling  challenge  is  to  express  damage  in  appropri¬ 
ate,  additive  units  for  situations  where  there  is  no  accepted  monetary 
surrogate  for  the  risk. 
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3.  NOTATION  AND  MOMENTS 

The  success  of  the  approximation  procedure  to  be  described  depends 
upon  the  assumption  that  most  of  the  total  risks,  (x^)  ,  have  a  high 
probability  of  being  zero;  this  can  occur  either  because  f°(0)  is  large 
or  because  n^  is  often  zero.  We  now  change  to  a  traditional  notation 
(see  JS)  which  emphasizes  the  distribution  of  risk  when  it  is  positive. 
Let 

(3.1)  Pr  {itj  -  0}  -  f°(0)  -  pt  -  1  -  q±  ,  (i  -  1,2 . N) 

(j  *  1,2,  •••» 

(3.2)  f±(x)  -  Pr  (x^  *  x  |  >  0}  -  f°(x)/q±  ,  (x  =  1,2 . R) 

and  define  the  first  two  moments  of  non-zero  risk  as: 

(3.3)  mt  -  E{xaj  |  x^j  >  0}  -  ECx^}/^  , 

(3.4)  v1  -  Vtx^  |  x^  >  0}  -  [Wx^J/q^  -  p^n^)2  • 

From  the  joint  counting  density,  we  get  the  marginal  densities: 

(i  -  1,2 . N) 

(3.5)  it  (n)  -  Pr  (n.  -  n) 

(j  -  1,2 . i±) 

and  the  first  two  moments: 


(3.6) 

(3.7) 


X^  ■  E{n^}  , 


Tik  ■  c<vV  • 


(i  ,k  =  1,2,  ...,  N) 


The  approximation  itself  is  based  upon  moment-matching  with  the 
first  two  moments  of  the  exact  density  (1.5),  which  we  now  find  in  a 


straightforward  manner.  First,  from  (1.3)  and  the  assumptions: 


E{xi  |  n^  ■  niE{x1j}  , 


(n/(i..}  (i  *  k) 

C{x  ;jL  |  n  ,n,  }  -  <  i 

1  K  1  K  (o  an) 


so  that,  unconditioning,  we  have: 


(3.8) 


E(ii>  -  E{ni)E{ilj}  , 


(E{n  }l/{x. . }  +  V{n.}[E{iL,}]2  (i  -  k) 

(3.9)  1  ij  1  ij 

(0  +  C(ni;nfc}E{xij  JECxj^^}  .  (i  i  k) 

Then,  using  (1.4)  and  notation  defined  above,  we  find  the  first  two 
moments  of  total  portfolio  risk  as: 


(3.10) 


and 


N 


(3.11)  l/{y} 


E{y>  -  I 

i-i  111 


N  /  2v  N  N 

■  Ji  >iMvi +  pi"i) +  j,  Yik’iVi”t 


The  (q^)  ,  (m^)  ,  and  (v^)  are  presumed  known  from  past  portfolio 

statistics  on  each  type  1  ,  and  the  (X^)  and  Cy^)  are  gotten 

from  modelling  assumptions  regarding  the  future  composition  of  the 

portfolio;  so,  we  shall  assume  that  these  moments  are  given  parameters. 

Of  course,  if  the  (n^)  are  statistically  independent,  the  last  term 
r  2  2 

become  l  "f^q^m^  •  *n  the  static  portfolio  model  in  JS,  the  composition 
was  fixed,  with  all  n^  •  1  ;  an  equivalent,  but  slightly  generalized 
model  can  be  gotten  from  the  above,  with  n^  -  n^  -  X^  ,  and  all  *  0  . 


'  «*  «v 


9 


.  THE  APPROXIMATING  RISK  COLLECTIVE  MODEL 


In  the  approximation,  we  replace  the  original  portfolio  by  a 
homogeneous  "risk  collective,"  that  is,  we  assume  that  y  is  approxi¬ 
mately: 


(4.1) 


a  =  o) 

e 


v  +  *  +  ...  +  w~  ,  (n^  >  0) 

e 


where  ng  is  the  random  number  of  equivalent  positive  claims  (w^)  , 
assumed  to  be  independent  of  each  other  and  ,  and  identically 
distributed,  according  to  prototypical  counting  and  individual  risk 
densities: 


(4.2)  ir(n)  *  Pr  {n  *  n}  (n  =  0,1,2,  ...)  ;  f(vv  =  Pr  {w  =  w}  (w  «  1,2 
leading  to  the  usual  compound  law  of  risk  theory  for  the  density  of  y  : 


(4.3) 


g(y)  ■  I  *(n)[f(y))' 
n-0 


As  mentioned  earlier,  the  rationale  behind  this  approximation  is  that, 
in  many  applications,  the  (x^)  are  zero  with  high  probability;  the 
(w^)  then  represent  just  the  positive  (x^)  .  (See  also  JS  and  Gerber 
(1979).) 

If  the  prototypical  moments  are: 


Chen  Che  momenCs  of  Che  random  sum  in  Che  approximaCing  model  will  be: 

(4.6)  E{y}  ■  Xm  , 

(4.7)  V{y}  ■  Xv  +  ym^  . 

For  a  good  approximacion,  Che  momenCs  (4.6),  (4.7)  musC  be  marched  as 
closely  as  possible  wich  Che  Crue  values  (3.10),  (3.11).  In  addicion, 
Che  forms  of  Che  ir(n)  and  f(w)  chosen  may  also  be  varied. 


5.  THE  ADELSON-PANJER  RECURSIVE  ALGORITHMS 

At  this  point,  we  should  stop  and  consider  whether  the  computa¬ 
tion  of  the  compound  law  (4.3)  can  be  effected  in  any  efficient  manner; 
otherwise,  it  is  not  much  improvement  over  (1.5).  A  traditional  approxi¬ 
mation  (for  the  static  portfolio  problem)  used  in  actuarial  circles  was 
to  make  x(n)  a  Poisson  law;  this  was  because  (further)  approximations 
to  the  compound  law  had  been  developed  in  the  early  risk  theory  litera¬ 
ture  (see,  e.g.,  Gerber  (1979)). 

However,  the  recent  extension  by  Panjer  (1981)  of  a  recursive  scheme 
of  Adelson  (1966)  now  provides  an  efficient  and  direct  way  to  compute 
(4.3).  Essentially,  if  f(w)  is  discrete  over  [1,2,  ...]  and  the 
counting  distribution  is  chosen  from  a  certain  (a3b)-family  for  which: 

(5.1)  ir(n)  -  ^a  +  ^)ir(n-l)  ,  (n  -  1,2,  ...) 
then  g(y)  can  be  calculated  recursively  via: 

g(0)  -  it (0)  -  (1  -  a)((a+b)/a) 

(5.2)  min (y , R) 

g(y)  -  1  (a  +  b  — )f (x)g(y  - x)  .  (y  -  1,2,  ...) 

x«l  '  y/ 


This  is  clearly  an  efficient  computational  procedure,  provided  the 
(a,b)-famlly  is  a  useful  one.  As  elaborated  upon  in  Sundt  and  Jewell 
(1981),  the  only  members  of  this  family,  apart  from  the  degenerate 
density,  are: 

n  -X 

(5.3a)  (Poisson)  ir(n)  -  —  (a  -  0  ;  b  -  X)  ;  (g(0)  -  e  D)  ; 


(5.3b)  (Binomial) 


(a  -  -p/(l-p)  ;  b  -  (-a)  (M  + 1) ) 


*(n)  "  wp  (1_P) 


(5.3c)  (Negative  Binomial)  ir(n)  *  (^a  ~  1jp°(l  - p)a  (a -p  ;  b  «p(o -1)) 

(Of  course,  the  familiar  Geometric  density  is  obtained  with  a  ■  1  in 
(5.3c).)  These  counting  distributions  are  useful,  since  they  are  often 
used  in  modelling  compound  risk  laws.  Furthermore,  since: 


(5.4) 


E{n  J 

e 


a  +  b 
1  -  a 


5  Y  =  Wne> 


a  +  b 
(1  -a ): 


we  get: 


(5.5) 


b  ,  x-CL.t  .1) 

Y 


and: 


(5.6) 


v  i 

Eii  )  X  1  ‘  s 
e 


The  importance  of  the  ratio  (5.6)  in  modelling  empirical  counting  processes 
is  well  known.  From  (5.3),  we  see  that  this  family  covers  a  wide  range 
of  such  ratios,  with  the  Binomial  giving  (yM)  <  1  and  the  Negative 
Binomial  (Pascal)  giving  (yA)  >  1  ;  the  Poisson  (y  “  1)  distribution 
is  the  dividing  line. 

Therefore,  for  computational  simplicity,  we  propose  to  use  the  (a,b)- 
family  to  model  the  counting  distribution  ir(n)  and  the  recursive  proce¬ 
dure  (5.2)  to  compute  the  approximate  density  (4.3).  Note  that,  if 
a  <  0  in  (5.5),  we  are  not  completely  free  in  our  choice  of  b  ,  since 
M  must  be  an  integer  in  the  Binomial  law  (5.3b);  however,  this  is  not 
usually  a  serious  limitation  (see  JS). 
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6.  THE  FINAL  APPROXIMATION 

Having  selected  n(n)  on  the  basis  of  computational  convenience, 
we  must  now  choose  the  prototypical  density,  f(w)  .  The  form  which 
will  give  the  best  approximation  in  all  cases  is  not  known.  However, 
the  most  natural  way,  consistent  with  the  interpretation  given  in 
Section  4,  is  to  weight  the  individual  densities  (3.2)  with  their 
average  chance  of  being  selected  and  giving  a  positive  ,  i.e., 

to  fix: 


(6.1) 


f  (w) 


I  Xiqif1(w)  l  X4fJ(w) 

I  Yj  l  Vj 


(w  -  1,2,  . . . ,  R) 


This  choice  is  consistent  with  JS  for  the  static  risk  portfolio  model, 
and  also  provides  the  greatest  simplification  to  the  formulae  below. 

Using  (3.3),  (3.4)  in  (6.1),  we  find  first  m  and  v  in  (4.5),  then 
substitute  into  (4.6),  (4.7)  to  find  the  first  two  moments  of  the  approxi¬ 
mating  model;  these  moments  are  then  equated  with  the  exact  results  (3.10), 
(3.11),  obtaining  finally  the  first  two  moments  of  the  prototypical  count¬ 
ing  density  in  terms  of  the  original  parameters : 


N 

(6.2)  A  -  E{n  )  -  l  A  q  ; 

i«l  1  1 


(6.3) 


Y  -  V(Se} 


N  f  /m  \2' 

'  i-i 


N 

+  l 

i-1 


N 

l 

k-1 


Yikqiqk 


m 


where  the  mean  prototypical  severity  is: 


l  Mi-i 
I  Yj 


(6.4) 


m  ■  E{w) 


and  the  severity  variance  is: 


To  summarize:  In  the  final  approximation,  we  would  first  calculate 
the  f^(x)  and  the  moments  of  Section  3  using  the  data,  then  compute 
f(w)  from  (6.1)  and  use  it  in  the  approximating  model  (4.3),  together 
with  one  of  the  ir(n)  of  Section  5,  with  (a,b)  selected  using  (5.5); 
the  approximate  density  is  computed  recursively  via  (5.2). 

In  the  static  portfolio  case  considered  in  JS,  all  y  are  iden¬ 
tically  zero,  so  that  y  <  A  ,  and  a  Binomial  counting  law  results. 

This  raises  the  integrality  problem  for  M  previously  mentioned,  and 
means  that  the  resulting  values  of  (a,b)  do  not  exactly  match  l/{y} 
in  the  original  and  approximating  models;  however,  the  resulting  error 
is  not  serious  in  the  example  analyzed  in  that  paper. 

In  contrast,  the  dynamic  portfolio  model  of  this  study  can  give 
y/A  >  1  ,  and  hence  Negative  Binomial  tt  (n)  ,  if  the  (y^)  are  large 
enough.  To  see  this,  consider  the  case  of  independent ,  but  still  random 
(n^)  .  (6.3)  then  becomes: 


(6.6) 


Thus,  we  see  that,  if  a  sufficient  number  of  (marginal)  counting  densi¬ 
ties  (3.5)  have  >  1  »  then  also  y/A  >  1  ,  a  most  reasonable 


result. 
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7.  THE  COMPOUND  MULTINOMIAL  COUNTING  DISTRIBUTION 

One  natural  way  In  which  the  number  of  risks  in  the  different  classes, 
n  -  (n^,^,  ...»  rijj)  ,  might  be  generated  in  a  predictive,  dynamic  model 
is  from  a  Multinomial  law,  with  given  total  number  of  risks,  n^  ,  and 
a  set  of  selection  probabilities ,  tt  -  ...,  irN)  ,  viz: 


With  fixed  n,j,  and  tt  ,  there  are  already  correlations  between  the 
counts  in  different  classes,  as: 


(7.2) 


(7.3) 


|  n^.:*}  -  ;  (i  -  1,2 . N) 


2 

Vt  -  W 
-  Vk°T  ' 


(i  -  k) 
(i  *  k) 


However,  to  give  more  modelling  flexibility,  we  now  permit  both  rtj.  and  tt 
to  be  a  random  scalar  and  random  vector,  respectively,  but  require  that 
they  be  independent  of  each  other,  for  simplicity.  This  "collective" 
model  dependency  gives  a  more  complex  covariance  structure. 

Define: 


N  __  N  N 

(7.4)  E^}  -  xT  -  J  ;  Wij)  -  yt  -  Yik  ; 

then,  unconditioning  (7.2),  (7.3),  we  obtain  the  moments  for  use  in  (6.2), 
(6.3): 

(7.5)  E(ni)  -  \±  -  XTE{iri)  ; 


<7. 6)  «S1;5k>  -  rlk  - 


XTE^i^  +  (yT -XT  +Xx)J/^1Ti^  +^YT  ~Xt^  t1 


(YT  “XT +Xl)C^i;^j^  +  (yt  -  XT) E{ tt± }E{ 7Tj  }  .  (i  i 


It  is  easy  to  show  that  these  satisfy  (7.4),  by  using  £  if  ■  1  . 

It  seems  to  the  author  that  practical  modelling  variations  might  fall 
into  one  of  two  extremes:  either  (1)  the  (tt^)  might  be  known  rather 
precisely,  and  forecasting  uncertainty  might  be  associated  with  the  total 
number  of  risks  or,  (2)  there  would  be  a  relatively  stable  number  of  risks 
but  prediction  uncertainty  would  remain  about  their  distribution  over  the 
different  risk  classification  types.  (For  the  casualty  claim  model,  only 
the  first  variation  would  probably  be  relevant.) 

An  interesting  special  case  of  the  compound  Multinomial  occurs  when 
the  (ir^)  are  fixed,  and  yT  “  XT  *  It:  then  follows  from  (7.5),  (7.6) 
that  A^  «  (i  -  1,2,  ...,  N)  ,  and  (n^n^i  4  k)  are  uncorrelated. 

This  then  simplifies  (6.2),  (6.3),  (6.6)  to  A  -  y  ,  that  is,  a  -  0  ,  b  - 
and  a  Poisson  counting  distribution  would  be  used  in  the  approximation 
of  Section  5!  One  obvious  way  in  which  this  could  happen  is  if  n^  were 
Poisson  with  parameter,  say  y  ;  it  is  then  well  known  that  the  (h^) 
must  be  statistically  mutually  Independent,  with  marginal  densities  that 
are  Poisson  with  parameters  (ir.,y)  . 


8.  AN  EXACT  RESULT 


There  is  one  case  in  which  the  proposed  procedure  gives  an  exact 
result.  Consider  a  risk  portfolio  of  fixed  size  N  ,  with  each  contract 
i  -  1,2,  N  having  an  Individual  claim  density  f°(x)  ,  with  parameters 

,  m^  ,  v^  ,  and  an  independent  claim  number  density  that  is  Poisson, 
with  parameter  .  This  is  the  basic  model  used  in  casualty  insurance. 

Following  the  procedure  in  Sections  5  and  6,  we  get  the  same  special 
results  described  in  the  previous  section,  namely,  y ^  ■  X^  ,  y  ■  0  (ii*k) 
and  X  -  y  -  I  *  *n  ot^er  words,  once  f(w)  is  determined  from  (6.1), 

the  recursive  algorithm  (5.2)  is  used  with  the  Poisson  density  (5.3a)  to 
find  the  approximate  g(y)  . 

However,  it  is  easy  to  show,  using  generating  functions,  that  the  exact 
form  (1.5)  reduces  to  a  compound  Poisson  law  with  parameter  X  ,  and  a 
severity  density  f(w)  .  Thus,  the  dynamic  portfolio  approximation  is,  in 
fact,  exaot  for  Independent  Poisson  claims.  This  is  true  even  if  p1  «  0 
for  all  1  i ! 

Unfortunately,  the  same  line  of  proof  shows  all  that  Independent 
Binomial  or  Negative  Binomial  claim  densities  (with  different  parameters 
for  each  i)  can  only  lead  to  an  approximation  of  the  true  g(y)  .  However, 
it  follows  from  Section  6  that  the  approximating  law  for  ng  would  be 
Binomial  or  Negative  Binomial,  respectively. 
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MODELLING  WITH  FIXED  AND  RANDOM  NUMBER  OF  COUNTS 


To  highlight  the  differences  between  the  model  and  procedure  of  this 
paper  and  the  static  portfolio  model  in  JS,  it  is  instructive  to  re¬ 
examine  how  the  independent  Poisson  casualty  claim  model  of  Section  8 
would  be  handled  according  to  the  JS  procedure.  We  use  primes  to 
designate  the  equivalent  parameters  of  this  paper,  in  terms  of  the  given 
model  parameters  y^  ,  q^  ,  m^  ,  v^  . 

First  of  all,  since  all  n^  =  1  in  the  JS  model,  we  would  have  to 
estimate  or  calculate  separately  the  N  individual  total  severity  den¬ 
sities  for  each  contract  risk,  x.  : 


(  )ne"Vi 

(9.1)  f°'(x)  *  g±(x)  *  l  — -  |f°(x)Jn*  * 


(This  could  be  done  by  N  applications  of  the  Adelson  algorithm,  or 
might  be  approximated  from  real  total  severity  data.) 

Then,  in  terms  of  the  parameters  of  this  paper,  we  would  get: 


(9.2) 


_yiqi  /qi\ 


Thus,  the  static  portfolio  approach  of  JS  would  use  the  Panjer  recursive 
algorithm  with: 

(9.3)  f'(w)  -  f°'(w )j/|l  qjj  »  (w  -  1,2,  ...) 


and  a  Binomial  counting  density  with  moments: 


vSO- .7  \v, 


9. 


10.  OTHER  VARIATIONS 


In  JS,  an  Improved  approximation  for  the  example  considered  was 
obtained  by  modifying  the  ir(0)  of  the  Binomial  (5. 3b)  to  enable  an 
exact  match  of  l/{y}  ,  together  with  an  integral  value  of  M  .  This 
modification  could  be  used  with  the  model  of  this  paper  whenever 
(yA)  <  1  ,  and  requires  only  a  trivial  change  in  the  recursive  algorithm. 

But  this  refinement  is  not  necessary  in  the  other  cases,  as  V{y}  is 
matched  exactly.  Of  course,  one  might  try  matching  other  moments  or 
values  of  the  exact  distribution  by  modifying  the  initial  values  of  the 
prototypical  counting  density  (see  the  discussion  in  JS). 

It  would  also  be  desirable,  particularly  in  Investment  applications, 
to  extend  the  range  of  permitted  (x^)  to  negative  values.  The  diffi¬ 
culty  then  is  that  the  relationship  (5.2)  is  no  longer  recursive,  and 
must  be  solved  by  other  means,  such  as  iterative  methods.  This  point 
Is  discussed  in  Sundt  and  Jewell  (1981),  where  possible  procedures  for 
the  Binomial  and  Poisson  cases  are  suggested;  exact  recursion  with  negative 
values  in  the  Negative  Binomial  case  (y/X)  >  1  does  not  seem  to  be  possible. 
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11.  COMPUTATIONAL  CONCLUSIONS  AND  ACKNOWLEDGEMENT 

The  limited  computations  carried  out  thus  far  indicate  that  the  same 
general  kinds  of  approximation  error  result  as  in  JS;  in  other  words,  the 
underlying  severity  density  should  not  be  too  "lumpy"  if  there  are  only 
a  few  risk  types.  Errors  also  seem  higher  in  strongly  correlated  cases, 
as  expected.  A  future  paper  will  explore  computational  results  in  more 
detail. 

The  author  would  like  to  thank  the  referee  who  found  several  errors 
in  the  original  formulae. 
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APPENDIX  A 


DEPENDENT  RISKS 

In  Section  1,  it  was  assumed  that  the  individual  risk  severities  (x^) 
were  statistically  independent  of  each  other  and  of  the  counts  (n^)  .  In 
this  appendix,  we  consider  the  modifications  necessary  if  the  risks  are 
exchangeable  random  variables  within  each  type  i  ,  but  still  independent 
of  the  counts.  As  is  well  known,  this  weak  dependency  is  equivalent  to 
assuming  that,  for  each  type  i  -  1,2,  ...,  N  ,  there  exists  a  random 
parameter,  0^  ,  such  that  the  individual  risks  are  independent  if  0^  =  0^ 
is  known,  and  depend  in  the  same  way  upon  0^  .  Thus,  the  basic  density 
(1.1)  is  replaced  by: 

(i  =  1,2,  ...,  N) 

(A.  1)  Pr  {x  =  x  |  0  }  -  f  (x  |  0.)  , 

J  1  (j  -  1,2,  ....  nj) 

giving  a  joint  density  within  type  i  ,  given  n^  *  n^  similar  risks,  of: 


(A.  2) 


Pr  <  n  X 

J-l 


ij  ij 


E  n 
j-i 


f?(x« 


and  a  common  marginal  density  for  any  risk  of  type  i  : 


(A. 3)  Pr  {x^  ■  x}  ■  Ef°(x  |  • 

(Expectations  in  the  above  are  over  the  random  values  of  Exchangeable 

random  variables  thus  have  the  property  that  they  have  the  same  marginal 
density  (and  self  moments) ,  their  arguments  may  be  permuted  in  any  fashion 
in  their  joint  density  (A. 2),  and  they  have  common  cross  moments. 


In  addition  to  the  dependency  between  different  types  introduced 

by  the  correlation  between  different  counts,  we  will  also  permit  the 

different  parameters  in  0  =  (0,,SO,  ....  0„)  to  be  statistically 

—  1  L  N 

dependent,  with  arbitrary  joint  d.f.  U(0)  .  In  short,  our  new  model 
substitutes  for  (1.5)  the  general  form: 


g(y) 

(A. 4) 


//  / 


dU(£)  l  l 


61  6  2 


nl  n2 


Intuitively,  we  can  think  of  0^  as  representing  exogenous  factors, 
such  as  the  economy,  weather,  political  factors,  etc.  that  influence  the 
random  outcome  of  all  risks  of  type  i  jointly.  This  type  of  "collective 
behaviour"  model  is  often  used  in  casualty  insurance,  where  it  is  recog¬ 
nized  that  all  risk  classification  schemes  are  imperfect,  and  that 
residual  correlations  still  exists  among  risks  of  a  given  type  due  to 
the  unexplained  inhomogeneity  still  present  within  the  class  i  .  Further, 
there  might  be  common  factors  between  the  different  classes,  which  would 


account  for  the  dependency  between  0^  and  0^  (i  /  k)  . 

Proceeding  in  a  manner  similar  to  Section  3,  we  define  the  positive 
risk  densities  f^(x  |  0^)  ,  the  probabilities  p^O^)  and  q^B^)  » 
and  the  first  two  moments,  m^ (0^)  and  v^ (0^)  ,  all  dependent  upon  the 
risk  parameter.  (3.8),  (3.9)  still  have  the  same  form,  except  that  they 


express  only  the  conditional  mean  total  risk,  E{x^  | 
covariance  of  total  risks  between  different  classes. 


0^}  ,  and  conditional 
|  e±;0k}  in 


terms  of  the  conditional  moments  of  individual  risk,  and  the  (non-0- 
dependent)  moments  (3.6),  (3.7)  of  the  counts. 


Now  all  that  remains  Is  to  uncondition  these  moments,  using  the 


relationships: 

N 

(A. 5)  E{y}  =  l  EE{x  |  8  }  , 

i-1  1  1 

N  N 

(A. 6)  l/{y}  *  I  I  ECCx^^  |  6^^)  +  C{E{x±  |  ;  E^  |  0k>}  . 

i»l  k=l 

(Innermost  operators  are  over  the  total  risks  (x^)  ;  outermost  operators 
are  over  the  risk  parameters  (0^).) 

We  define  the  unconditional  versions  of  q^(8)  ,  1^(0)  ,  v^(0)  as: 

(A. 7)  q±  =  E{qi(0)}  ;  m±  -  Etm^Oj)}  ;  v±  -  E{vi(§i)}  . 

By  the  theorem  of  conditional  expectation,  q^“  Pr (x^  >0}  is  the  same 
as  in  (3.1).  However,  as  the  referee  reminds  us,  m^  and  v^  are  not 
the  same  as  and  v±  in  (3.3),  (3.4)  unless  the  variation  due  to  0i 

vanishes;  hence,  the  different  notation.  In  fact,  in  the  current  notation 
we  see  that: 

(A. 8)  I  >  0}  -  mi  +  C{qi(§i)  ;  mi(01)}/qi  . 

In  addition  to  correlations,  we  shall  also  need  higher-order  cross- 
moments,  so  we  define: 

(A. 9)  Qi(0i)  -  qi(0i)  -  q±  ;  M1(0±)  -  m^)  -  ;  V±(01)  -  v±(01)  - 

and  use  notation  like: 
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■  EtVVVV>  -  ; 


(A. 10) 


QiMi  ’  «VVVV>  ! 


TO  ’  RVVWW1  > 


and  so  forth. 


In  place  of  (3.10),  we  have: 


(A. 11) 


E(y)  -  l  + 1 1  . 


and.  In  place  of  (3.11),  we  obtain: 


W)  -  I  +  Pi^,)2)  +  l  l  ^kWA 

l  -  AM) 

*  I  »i[vT  +  zijCPi  -  q^QX  +  <P±  -  Wi 
-  2S1qJm1  -  q^] 


(A. 12)  + 


+  l  J  2Yikqimi\Mk 

1  k 


+  l  £  <Yik +  W  [qiqkMA  +  ;iTO  +  !TOA 


2qiWk  +  fiiTO 


wa] 


[{  xlQlMl]: 


V-'.  - .  .  • 


The  term  in  braces  in  (A. 11)  gives  a  correction  term  to  the  calculation 
of  X  in  (6.2)  (with,  of  course,  m^  and  replaced  by  m^  and  v^); 

similarly,  the  terms  in  braces  in  (A. 11)  and  (A. 12)  give  two  correction 
terms  to  the  calculation  of  y  in  (6.3). 

In  many  applications,  these  corrections  simplify  because  either  the 
probability  of  a  claim  or  the  moments  are  independent  of  6^  .  For  instance, 
in  life  insurance,  m^  -  m^  and  v^  ■  v^  are  the  moments  of  the  face  value 
of  policies  of  type  i  ,  which  do  not  usually  change  with  exogenous  condi¬ 
tions,  while  the  expiration  probability,  q^(0)  ,  would  probably  vary  with 
external  effects;  this  would  eliminate  all  terms  in  (A. 11),  (A. 12)  with 

,  or  !  Conversely,  in  casualty  Insurance,  the  probability  of 
a  claim,  ,  might  be  relatively  fixed  several  years  in  a  row,  but  the 
severity  moments,  *1^(0^)  »  night  be  relatively  uncertain  in 

view  of  Inflation,  etc.;  in  this  case, all  terms  in  (A. 11),  (A. 12)  involv¬ 
ing  and  can  be  eliminated! 

A  more  complex  model  can  also  be  developed  by  permitting  the  (n^) 
to  depend  upon  ;  details  are  left  to  the  reader. 
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Y  “  Wne) 

\  ni 

it.  should  be  it. 
i  i 


page  17,  line  -5 


"all"  should  be  "also" 
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Correction 


p.  9,  line  -5  and  -6 


p.  9,  eq.  (4.4) 


p.  11,  eq.  (5.2) 


p.  11 ,  eq.  (5.3a) 


p.  12,  eq.  (5.3b) 


p.  12,  first  paragraph 


p.  13,  first  paragraph 


p.  15,  in  (7.1) 


p.  17,  line  -5 


(x.)  should  be  (x. .) 

B  I  J 


Y  =  V{ne} 


g(0)  =  n(0)  =  (1  -  a)(a+b)/a  (a  t  0) 


U  =  0) 


delete:  (g(0)  =  e~b) 


b  =  -a(M  +  1) 


delete:  (Of  course,  the  familiar  ...  (5.3c).) 


However,  a  natural  way,  consistent  with  the  inter¬ 
pretation  given  in  Section  4,  is  to  weight  the 
individual  densities  (3.2)  with  weights  proportional 
to  the  expected  number  of  risks  with  positive  outcome 
in  the  corresponding  class,  i.e.,  to  fix: 

nk  ni 

ir.j  should  be 


'all"  should  be  "also" 


:*r it*  t 


p.  18,  eq.  (9.2) 


vi  ‘  (^f)  Vi  *  <mt  -  • 


